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Abstract: 

We consider the integrable spin chain model - the noncompact SX(2, R) spin magnet. The 
spin operators are realized as the generators of the unitary principal series representation of 
the SL(2,M) group. In an explicit form, we construct IZ— matrix, the Baxter Q— operator and 
the transition kernel to the representation of the Separated Variables (SoV). The expressions 
for the energy and quasimomentum of the eigenstates in terms of the Baxter Q— operator are 
derived. The analytic properties of the eigenvalues of the Baxter operator as a function of the 
spectral parameter are established. Applying the diagrammatic approach, we calculate Sklyanin's 
integration measure in the separated variables and obtain the solution to the spectral problem 
for the model in terms of the eigenvalues of the Q— operator. We show that the transition kernel 
to the SoV representation is factorized into a product of certain operators each depending on a 
single separated variable. 
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1. Introduction 

Solutions of many integrable models can be obtained within the Algebraic Bethe Ansatz (ABA) 
method |1J. However, in some cases, this method is not applicable. This happens, in particular, 
when the Hilbert space of the model does not possess the lowest weight vector. In this case the 
solution of the model should rely on the methods of the Baxter Q— operator [2j and the Sepa- 
rated Variables (SoV) jHlH], which present an alternative to the ABA approach. The advantage 
of these two methods is that their applicability is not subjected to the restrictions of the ABA 
method, while the disadvantage is that the Baxter Q— operator and representation of the Sep- 
arated Variables are known only for a limited class of models. The most well known model of 
this class is the quantum Toda chain. The method of SoV was developed by Sklyanin ,3J and 
worked out for the Toda chain by Kharchev and Lebedev [Hj. The Baxter Q— operator for the 
Toda chain was constructed by Pasquier and Gaudin 0. 

Recently the Baxter Q— operator and representation of the Separated Variables have been 
constructed for a number of models El El EE UH EH El CHI 13 EEj- Among these models 
there are so-called noncompact spin magnets. They are the cousins of the conventional spin 
magnets like the famous XXX Heisenberg spin magnet. The difference between these models is 
that the spin generators act in different Hilbert spaces. In the case of the compact magnets it 
is the finite dimensional space of certain representations of the SU (2) group. The Hilbert space 
of the noncompact magnets is taken to be the vector space of the unitary representations of the 
noncompact group SX(2,R) or SX(2,C), which are necessarily infinite-dimensional. 

Compact SU(2) spin magnets are of special interest in statistical physics (see e.g. Ref. fTKj). 
The recent interest in the study of the noncompact magnets originates from high energy physics 
(see review Refs. [T71 CHI)- The studies of the scale dependence of scattering amplitudes in 
Quantum Chromodynamics revealed the surprising fact that the equations governing the scale 
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dependence of many physically important amplitudes are integrable. Namely, it has been shown 
that the Hamiltonian governing the Regge-behaviour of the hadronic amplitudes is equivalent to 
the Hamiltonian of the noncompact SL(2, C) spin magnet 0120] • Further, it was observed in [2T] 
that the scale dependence of certain partonic distributions in QCD is governed by Hamiltonians 
which are equivalent to those of the spin chains, both closed and open, with SX(2,R) symmetry 
group. The Hilbert space of the latter is given by the tensor product of the vector spaces of 
the discrete series representations of the SL(2, R) group and possesses the lowest weight vector. 
Therefore the solution of these models can be obtained by the ABA method (see Refs. (22] 1221 EI] 
for the analysis of these particular models). Contrary, the solution of the SL(2, C) spin magnets 
relies entirely on the method of the Baxter Q— operator and Separation of Variables jTTj. 

A general method of constructing the Baxter Q— operator is not developed yet and the latter 
is known for a limited class of models only. For the noncompact SL(2, R) spin chain the Baxter 
Q— operator was first constructed in Ref. [Zj. The approach used in |Zj is based on the Pasquier 
and Gaudin method |6j and suggests an effective way to resolve the defining equations and obtain 
the integral kernel of the Q— operator in an explicit form. Further, in the case of the noncompact 
spin magnets it appears to be quite helpful to interpret the integral kernels of the operators in 
question (Q— operator, transfer matrix, etc.) as Feynman diagrams of a certain type [TTj I12j. 
Then the analysis of the properties of the models is drastically simplified and in many cases can 
be fulfilled diagrammatically. Moreover, following this approach one can construct not only the 
Baxter operator that allows to determine the energies of the eigenstates, but also the transition 
operator to the SoV representation [TTJ [T3J that gives the explicit representation for the 
eigenfunctions of the model. Surprisingly, such a reformulation appears especially effective for 
the models where the ABA method does not work, e.g. SL(2, C) spin magnet [TTj . 

In the present paper we consider the noncompact spin magnet with the Hilbert space given by 
the tensor product of the vector spaces of the unitary principal series representation of SX(2,R) 
group. Such models have not been considered so far and differ drastically in their properties from 
the SL(2, R) magnet considered in the literature. A particular choice of the Hilbert space gives 
this model some similarity with the SL(2, C) spin magnets. This similarity also appears in the 
method of the analysis used. The solution in both cases relies on the Baxter Q— operator and 
SoV methods. However, this model possesses some specific properties which make it different 
from both SL(2, C) and SX(2,R) (discrete series) spin magnets. The two-particle Hamiltonian 
has both the discrete and continuous spectrum, the energy of the corresponding eigenstate is not 
completely fixed by its two-particle spin and deviates from the usual ip— function form. We found 
that such fundamental object as 1Z— operator, Baxter Q— operator are doubled, i.e. there are two 
independent solutions of the Yang-Baxter relation for 1Z— operator and two Baxter Q— operators. 
The latter satisfy the Baxter equation which entangles them. Although it is possible to introduce 
a linear combination of the Baxter operators such that each satisfies the separate equation, they 
remain entangled through the quantization conditions. 

The paper is organized as follows: In the Sect. ED we remind some facts about SL(2, R) group 
and introduce the notations. In the Sect. |3Dwe construct the SX(2,R) invariant solution of the 
Yang-Baxter relation and discuss its properties. The transfer matrices and Hamiltonian are built 
in the Sect. 13.2.1 In the Sect. U]we construct the Baxter Q— operator and study its properties. 
The representation of the Separated Variables is constructed in the Sect ED The eigenvalues of 
the Baxter Q— operator can be found in analytic form for the two-site chain, which is considered 
in the Sect.|6D The concluding remarks are given in the Sect. [?D and the Appendix contains some 
technical details. 
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2. Preliminaries 



The aim of this section is to remind some basic facts about the representations of the group 
SL(2, R) and introduce the necessary notations. It is well known that all unitary representations 
of the group of the real unimodular matrices SL(2, R) can be organized into three series, the 
discrete ones, the principal and supplementary continuous series |2Hj . The latter will not appear 
in our analysis and therefore is not considered here. 

The unitary representation of the principal continuous series is determined by two numbers, 
a real p and a discrete e, which takes only two values and 1/2. It is convenient to denote the 
pair of numbers (a, e), where a is complex and e is or 1/2, by the bold letter ct. We define the 
sum of two numbers Qi and ct2 as follows 

Oil + «2 = «3 = («i + «2,e 3 ) • (2.1) 

where €3 = 0, if e% + €2 is integer, and €3 = 1/2 otherwise. In what follows writing the sum 
(ei + £2) we shall always imply such addition rule. Next, the zero element is defined as = (0, 0). 
By — ct we shall denote the element inverse to ct, ct + (—ct) = 0, — ct = (— a,e). If the first 
element, a, of the number ct is real, we shall call the number ct real as well. The usefulness of 
such notations will be clear later. 

Thus the unitary representation of the principal continuous series is labelled by the real p. It 
can be realized by the unitary operators T p (g) acting on the Hilbert space L 2 (M) |2H] 



= wrw^y^Td)' (2 - 2) 

where g _1 = ^ ^ j and the sign factor is cr e (x) = [sign(a;)] 2<E . All representations T p (except of 

the representation T^ ' 1 / 2 ) ) are irreducible. The two representations T p and T 9 ' are not equivalent 
provided that p 7^ —p'- The unitary operator Ai p which intertwines the representations T p and 
rp-p (t~ p A4 p = M. P T P ) is defined uniquely up to a phase factor 



[M p mx) = {^A^))- 1 fjyHv) |^| 2( f-l) . ( 2 - 3 ) 

where the parameter s (conformal spin) is determined as s = 1/2 + ip and 7 = (2 — 2s, e). The 
function A(a) is given by the following expression 

M a ) = ^ (f) 2e a(e + a/2), (2.4) 

where a(x) = T(l/2 — x)/T(x) . This function, which appears naturally in the course of the 
calculations, satisfies the following conditions 

A(a)A(l-a) = (-1) 2 % A(cx)* = (-l) 2 M(a*) , (2.5) 



where ct* = (a*, e). The normalization factor chosen in ()2.3)1 ensures that the intertwining 
operator satisfies the following equation A^J, = M.- p 1 . 

1 Henceforth we shall not display explicitly the label p of the operator M. since it completely fixed by the 
transformation properties of the function the operator M. acts on. 
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The three generators of the s£(2) algebra corresponding to the realization of the representation 
T p (EH have the form 



S- = —d x , S + = x 2 d x + 2sx, Sq = xd x + s . (2.6) 

These operators are antihermitean and obey the standard s£(2) commutation relations 

[S ,S ± ]=±S ± , [S + ,S^] = 2S . (2.7) 

The tensor product of two representations of the principal unitary series can be decomposed 
into the direct integral of the representations of the same type and the direct sum of the repre- 
sentations of the discrete series, V^, |26] 

dpT p + ( V t® V h)- (2-8) 

h=l+(e 1 +e 2 )/2 

The representation of the continuous series, T p , enters into the direct integral with the multiplicity 
two. The operators separating the irreducible components and the other necessary details can 
be found in the Appendix EJ 

We conclude this section with the following remark. Let us divide all tensor products T P1 ®T P2 
in two subsets depending on the value, or 1/2, of the sum e\ + e 2 . Then it is easy to see that 
all representations inside each group are unitary equivalent to each other. In particular, the 
representation of the first group, ei + e 2 = 0, are equivalent to T^ 0,0 - 1 ® T^°'°\ and those in the 
second group, e\ + e 2 = 1/2, are equivalent to T^ ' ) ® T^ ' 1 ^ 2 . To show this let us consider the 
operators V(ot), 

[V(ct)^)(x 1 ,x 2 ) = ° e{Xl _~ X *l ttfosa), (2-9) 

1^1 1 

and U(a), 

U(a) = (M <g> I) V{a) (M <g> I) , (2.10) 

where o: = (a, e) and the operator A4 defined in (12. 3j) intertwines the representations T Pl and 
T~ Pl . It is obvious that for real ct, the operators V(a) and U(ct) are unitary and that they 
intertwine the representation T Pl ® T p2 with T Pl+a ® T P2+a and with T Pl ~ a ® T P2+a , respec- 
tively. Therefore the combination, U(a) V(f3), with suitably chosen a,f3 intertwines any two 
representations inside each group. We notice also that for a real ct, 

V{ a y = V(-a), U( a y = U(-a) . (2.11) 



3. 1Z— operator 

In this section we shall construct the solution of the Yang-Baxter relation 

n 12 (u)n 13 (u + v)n 23 (v) = n 23 (v)n 13 (u + v)n 12 (u). (3.1) 

The operator 1Z ik in the above equation acts on the tensor product of the spaces Vi ®V k and 
depends on the spectral parameter u = (u,e), i.e. lZ(u) = lZ(u,e). Each space V$ is equivalent 

2 The representation T^ 0,1 / 2 ) is reducible and equivalent to the T>^, 2 
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Figure 1: The diagrammatic representation of the chain relation, Eq. ()3.9|) . The black dot denotes 
the integration vertex and 7 = a + (3 — 1. 
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Figure 2: The diagrammatic representation of the star-triangle relation (Eq. (j3.10|0 . The primed 
indices are defined as x' = 1 — x . 



to the Hilbert space L 2 (M) and carries the representation of the principal continuous series, T Pi , 
of the SX(2,R) group. Let us notice that Eq. (j3.1|) contains two equations involving operators 
R(u) = 1Z(u, 0) and R(u) = 1Z(u, 1/2). Indeed, having put u = (u, 0) and v = (v, 0) in Eq. (j3.ll) 
and taking into account that u + v = (u + v,0), one gets the Yang-Baxter equation for the 
operators R. At the same time, choosing u = (u, 1/2) and v = (v, 1/2) one finds that the 
argument of the operator 7li3 remains the same, u + v — (u + v, 0), and therefore, in this case, 
Eq. (|3.1|) involves two R operators and one R(u + v). At a given choice of the spectral parameter 
these two equations are contained in one equation (|3.1|) . It will be seen later that such, slightly 
unusual, choice of the spectral parameter is quite natural for this model. 
We shall look for the solution of ()3.1|) in the form of the integral operator 



dy x dy 2 R u (x 1 ,x 2 \yi,y2) * (2/1, 2/2) 



(3.2) 



and impose the additional restrictions of the SL(2, IR) invariance. The relation (|3.1|) leads to cer- 
tain equations on the integral kernel. We shall propose an ansatz for the kernel R u (xi, x 2 \yi, y 2 ), 
which is motivated by the form of the kernel of the 1Z operator for the SL(2, C) magnet [TT] . an d 
show that it satisfies Eq. (|3.1|) . To write down the kernel in a compact form let us define the 
function (propagator) D a (x), 

aAr) (3.3) 



x 



x 



which depends on real x and has the index a. 
the product of four propagators 



(a, e). The kernel of the 1Z— matrix is given by 



Ru(xx,x 2 \yi, y 2 ) 

3re r PiP2 (u) is th 
parameter. The requirement of the SL(2 



r PlP2 (u) D a4 (y 2 - xi) D a3 (xi - x 2 ) D a2 (x 2 - y x ) £>«i(?/i 



2/2) , (3.4) 



where r PiP2 (u) is the normalization coefficient to be defined later and u 



[u, e) is the spectral 
invariance of the kernel imposes the following 
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Figure 3: Diagrammatical representation of the kernel of the 7£-operator, Eq. ()3.4|) . 



constraints on the indices oti 

«i + «2 = (2 - 2a a , ei) , a 3 + at 4 = (2s lt e x ) , 

Oil + 0:4 = (2 - 2s 2 , e 2 ) , ck 2 + e* 3 = (2s 2 ,e 2 ) . (3.5) 

Therefore it is sufficient to fix only one index to restore all others. Again, using similarity with 
the SL(2,C) magnet (see Ref. [TT]) we put 

014 = (1 + ipi — ip 2 — iu, e) (3.6) 

that gives rises to the following values of the other indices 

oi\ = (-ipi - ip 2 + iu,e 2 + e) , 

ol 2 = (l — ipi+ip2 — iu,e 1 + e 2 + e), (3.7) 
«3 = (ipi + 1P2 + iu, ei + e). 

The proof of the Yang-Baxter relation is based on the integral identity for propagators - the 
so-called star-triangle relation. It is well known in the perturbative Quantum Field Theory and 
is widely used for the evaluation of the multi-loop Feynman diagrams [2Zj • As will be seen later, 
it is quite natural to represent the integral kernels of the operators under consideration in the 
form of Feynman diagrams. It allows one to carry out an analysis of the properties of the model 
diagrammatically, which in many cases results in considerable simplifications. Therefore, below 
we collect some elements of the "diagram technique" we use in our analysis, 

• Fourier transform 

J dxe ipx V a {x) = v / ^2 1 - a A(a)P 1 _ a (p), (3.8) 

where 1 = (1, 0) and the function A(a) (a = (a, e)) is defined in Eq. (|2.4J1 . 

• Chain relation 

J dyV a (x - y)V p {y - z) = v ^^^ V y (x - z), (3.9) 
where 7 = cc + (3 — 1. 
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Figure 4: The diagrammatic proof of the Yang-Baxter relation, Eq. 1)3.1)) . 

• Star-triangle relation 

J dwV a {x — w)Vp(y — w)V~f(z — w) = (3.10) 

= y/wA(a, (3, 7) Di_ Q (2 - y)T> 1 _ f3 (x - z)Vx^{y - x). 

Here A(ct,(3,j) = A(a)A((3)A('y) and indices a,/3, 7 satisfy the uniqueness condition 
a + /3 + 7 = (2,0). 

The diagrammatic representation of the above identities is given in Figs. ^ and 121 There the 
arrow directed from point y to x and labelled by the index a denotes the propagator T) a {x — y) 
and the black dot is used for the integration vertex. 

At last, we give the following representation for the delta function 

5(x) = lim^li^. (3.11) 

which can be obtained from Eq. (J3.8)) . 

The proof of the Yang-Baxter relation for the TZ— operator defined in Eq. (J3.4)) can be carried 
out diagrammatically. To this end let us notice that graphically the 1Z operator is represented 
by the box diagram (see Fig. EJ). Then the lhs and rhs of the Yang-Baxter relation ()3.1)) are 
represented by the diagrams shown in the the Fig. left up and lower diagrams, respectively. 
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Figure 5: The diagrammatic representation of the rhs of Eqs. (j3.14|) . 

To prove the Yang-Baxter equation one has to show that these diagrams are equal. The proof 
repeats the one given in Ref. ^T] and is based on the use of the star-triangle relation (j3.10|) . 
Indeed, an examination of the indices of the central triangles in those diagrams shows that they 
satisfy the uniqueness condition and, therefore, can be replaced by the star diagrams ()3.10p (see 
second up and lower diagrams in Fig. 0J). The new triple vertices in each diagram appear to be 
unique as well. In the next step one replaces the corresponding star subdiagrams by triangles. 
The resulting diagram, in both cases, is the hexagon diagram. Restoring all factors which arise 
in the course of above transformations, see Eq. ()3.10|) . one finds that the final diagrams coincide. 

One can also show that 7Z— operator satisfies the Yang-Baxter relation involving two Lax 
operators 

L x {u)L 2 {v + u)K 12 {u) = 1Z X2 {u) L 2 {v) L^u) , (3.12) 

where, as usual, Li{u) = L{u) ®I, L 2 (u) = I®L(u). The Lax operator is defined in the standard 
way 

L{u) = u + if |° _f~ o ^ (3.13) 

and is used for the construction of the auxiliary transfer matrix. 

Let us note that the above Yang-Baxter relation holds for the both values of the e, which 
enters the spectral parameter of the 1Z— operator, i.e. Eq. ()3.12|) holds both for 7Zi 2 (u,0) and 
K 12 {u,l/2). 

3.1. Properties of 1Z— operator 

Since the tensor products T Pl <g> T P2 fall into two classes of equivalence, it is natural to expect 
that the 1Z operators should be unitary equivalent as well. Indeed, it follows directly from the 
definition of the 1Z operator ()3.4|) that 7Z Pl + a , P2 +a ~ V( a )7Z PlP2 ^( — oc). In the general case 
the following relation holds: 

1Z PlP2 {u) = U{a)V(J3) 1Z m , {0 , ei+e2) (u, 0) V(-(3) U{-cx) , (3.14) 

where ex = ((p 2 — pi)/2,e) and j3 = ((p 2 + pi)/2,ei + e). Here some comments are in order. 
The diagram corresponding to the integral kernel of the rhs of Eq. (I3.14J) is drawn in Fig. |S| 
The central box corresponds to the kernel of the operator V(/3)7Z^ ^^o :tl+e2 )(u,0)V(—f3), and 
the star subdiagrams from the left and right sides represent the kernel of the U operators. The 
straightforward application of the star-triangle relation fl3.10|) in a way indicated in Fig. El - one 
applies the star-triangle relation to the vertices (a) and then to the vertices (b) - turns the 
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diagram shown in Fig. El into the box diagram with indices given by 1)3.6)1 . 1)3.7)1 . Therefore, up 
to some numerical factor, the resulting diagram coincides with those for the TZ— operator. 

Let us remind that the normalization factor, r PlP (it), in Eq. ()3.4j) is so far arbitrary. Since an 
arbitrary TZ— operator is related by Eq. (J3.14)) to the operator 72(o,o),(o,o)( M ) 0) or 72.(o,o), (0,1/2) {u, 0), 
we fix the normalization of the latter, and for all other operators take Eq. (j3.14j) for a definition. 
Namely, we fix the normalization factor in Eq. ()3.4)1 for the two selected operators as follows 

r(o,o),(o, ea )(«,0) = 2- 2iU — -A{iu,0)A{iu,e 2 ), (3.15) 

where the function A(ct) is introduced in (|2.4jl . Then one finds that the definition of the 
TZ— operator via Eq. 1)3.14)1 is equivalent to the following choice of the normalization factor in 
Eq. (|331) 

r PlP2 (u) = ^e l ^ ei+e2) 2- 2iu A(ip 2 -ip 1 + iu,e)A(ip 1 -ip 2 + iu,e + e 1 + e 2 ). (3.16) 

We remind that the sum (ej + e 2 ) is defined by Eq. 1)2.1)1 . Such a choice of normalization ensures 
that the TZ— operator is unitary for the real u, 

TZi2{u)TZ[ 2 {u) =1. (3.17) 

Using Eq. 1)3.11)1 it can be easily checked that for p 1 = p 2 

TZ 12 (0) = (-l) 2ei Pi 2 , (3.18) 

where P 12 is the permutation operator, [P 12 ty](xi, x 2 ) = ty(x 2 ,xi). 

Moreover, taking into account that V(—f3)U(—a) = (U (a)V ((3))\ (see Eq. (12.11)1 ). one 
concludes that Eq. ()3.14)1 establishes the unitary equivalence between the arbitrary TZ— operator 
and the operator 72.(0,0), (0,0) ( n ) ( or ^(0,0), (0,1/2) ( u ))- Thus in order to determine the spectrum of an 
arbitrary operator TZ PlP2 it is sufficient to calculate the spectrum of the TZ operator for the special 
values p 1 = (0, 0) and p 2 = (0, e 2 ). The eigenfunctions of the TZ operator are fixed by the SL{2, R) 
invariance of the latter. To find them one has to construct the operators separating the irreducible 
components in the tensor product of two representations of the continuous series, Eq. ()2.8jl . 
The tensor product decomposition ()2.8)1 contains the representations both of the discrete and 
continuous series. The operators (U. p ' p (xi,x 2 ,w)) separating the representations of the discrete 
series, T>^, are given in Eq. (JA.9)) 3 . Next, since the representations of the continuous series enter 
the tensor product decomposition Eq. ()2.8jl with double multiplicity, there is arbitrariness in the 
choice of the projectors. We defined the projectors to the continuous series, (xi, x 2 , y), 4 

(see Eq. ()A.4)l ) in such way that they possess definite parity with respect to the interchange of 
the arguments x\ and x 2 . A direct check shows that namely such combinations diagonalizes the 
TZ— operator. Introducing the notations for eigenvalues of the TZ— operator on the corresponding 
eigenfunctions 

TZ PlP2 {u,e) (llg^y = J^(u,e)(n£ ft )*, (3.19) 
(n^) 4 = J$(u,e)(n£^)*, (3.20) 

one obtains after some calculations 

3 We remind that h takes integer or half-integer values. 
4 The parameter e = 0, 1/2 marks different projectors. 
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Figure 6: The diagrammatic representation of the kernel of the transfer matrix, Eq. ([3.25)1 . The 
indices cxi of the lines on the fc-th box are given by Eqs. ()3.6j) and ()3.7j) with p l = p and 

P2 = Pk- 



ei = £2 



R ps (u,e) = - \sm(inu) + (-l) 2e+2£ sm(ns)}T(l~s-iu)T(s-iu), (3.21) 

7T 



(3.22) 



ei 7^ e 2 



R pE (u,e) = - \cos(i7ru) + (-l) 2e+2£ cos(7fs)]T(l-s-iu)T(s-iu), (3.23) 

' 7T 



P±/„ c ^ _ C ^fe-l/2 r (^ m ) 

r(/i + nt) 



(3.24) 



where s — 1/2 + ip. 
3.2. Transfer matrices and Hamiltonians 

Having obtained the solutions of the Yang-Baxter relation ()3.1jl . one can construct an infinite set 
of the commuting 5'L(2,R) invariant operators (transfer matrices) 

T Po (u) = tT PQ (n PoPl (u)...TZ PoPN (u)) , (3.25) 

where the trace is taken over the auxiliary space V Po . They depend on the spectral parameter u 
and the spin of the auxiliary space, p . The diagrammatic representation of the transfer matrix 
is shown in Fig. El Invoking the standard arguments |2Hj, one finds that the transfer matrices 
form the family of the mutually commuting operators 

[T p »,T p[) W] = 0. (3.26) 

In the case of the homogeneous chain, (when all quantum spaces carry the same representation 
of the SX(2,M) group), i.e. when p x — p 2 — . . . — p N = p q , it is possible to construct the 
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operator (Hamiltonian) which has two-particle structure. Indeed, choosing p = p q and taking 
into account the property (j3.18j) . one obtains 



n 



N 



H, , + ■ • • + H 



N~l,N 



+ H 



NA 



(3.27) 



u=(0,0) 

where the two-particle Hamiltonian H^^+i is given by 

d 



H 



12 



u=0 



(3.28) 



Introducing the notation E £ (p) for the eigenvalues of the pairwise Hamiltonian corresponding to 



the eigenfunctions IKl^ 5 and E ± (h) for those corresponding to the eigenfunctions n 



derives from Eqs. (|3.21|) the following expressions for the energi 



one 



ICS 



£ £ (p) 



V> ( - + ip 
2ip{h) . 



+ '*<2 



1 P 



\2e 



7T 



cosh 7rp ' 



(3.29) 
(3.30) 



Herein ip(x), as usual, denotes the logarithmic derivative of the Euler T— function. Thus one 
sees that the energies corresponding to the discrete levels are double- degenerate. The separation 
between the continuous branches being large in region of small p, vanishes rapidly for large p. 
The corresponding dispersion curves are shown in Fig. [71 Note also, that although the eigenvalues 
of the two-particle Hamiltonian do not depend on the p q , this dependence reveals itself in the 
eigenfunctions. For the chains with the number of sites N > 2 the energies will explicitly depend 
on the values of p q . 

Let us consider now the simplest among the transfer matrix — the auxiliary transfer matrix. 
It is given by the trace of the product of the Lax operators (|3.13j) 



t N (u) = tr (Li(u) . . .L N (u)) . 



(3.31) 



The auxiliary transfer matrix, ijv(it), is the polynomial of degree N in u with operator valued 
coefficients 



A' 



t N (u) = 2u N + J2 



u N - k q k 



k=2 



The operators qj~ (integrals of motion) are differential operators in Xk- In particular, 



S 2 - N(p 2 q + 1/4) 



(3.32) 



(3.33) 



where S = Si + . . . + S^- is the operator of the total spin. The possible values of the total 
Casimir operator, S 2 , are restricted by the SL(2, M) invariance, namely, S 2 = h(h — 1) for the 
representations of the discrete series and S 2 = — 1/4 — p 2 for the representations of the continuous 
series. 

Further, by virtue of Eq. ()3.12|) one derives that 



S, t N (u) 



[t N (u),t N (v)] = [t N (u),T Po (v))] = [t N (u),H N ] 



0. 



(3.34) 



3 Since p 1 = p 2 , one gets that p — (p, 0). 
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Figure 7: The dispersion curves for the two-particle Hamiltonian. The curves for E^(p) and 
E^ l ' 2 \p) are shown by the solid line and the dot-dashed line, respectively. The discrete eigenval- 
ues E ± {h) (h = 1, 2, ..) are depicted by the black diamonds. 

It follows from ()3.34|) that the auxiliary transfer matrix and Hamiltonian share a common set of 
eigenfunctions. Thus the eigenfunctions of the Hamiltonian can be labelled by the eigenvalues 
of the integrals of motion q 2 , . . . , qN and one of the SL(2, R) generators, which is usually chosen 
to be S-. However, in the case under consideration, these quantum numbers do not fix the 
eigenfunction completely. Indeed, as seen already for the two-site spin chain, the eigenfunctions 
of the Hamiltonian, corresponding to E°(p) and E 1 ^ 2 (p) (see Eq. (13.21)1 ) have the same quantum 
numbers. Therefore, the spectrum of the transfer matrix is double degenerate, and fixing of 
the total momentum, iS_ = p, does not remove this degeneracy. Nevertheless, the number 
of eigenfunctions with identical integrals of motion, q^, is finite. One can always resolve the 
degeneracy on this finite subspace by fixing an additional quantum number, for example the 
value of the quasimomentum 6. We remind that e is the eigenvalue of the operator of the cyclic 
permutation 

V^(xi, x 2 ,..., x N ) = $(x 2 , ...,x N , xi) , (3.35) 

which commutes with both auxiliary transfer matrix and Hamiltonian, [P, tv(w)] = [P, 'Hn] = 0. 
Thus to uniquely determine the eigenstate, one should the specify the total momentum, p, and 
the following set of the quantum numbers, q = {9, q 2 , ■ ■ ■ , (/at}. 

4. Baxter Q— operator 

The solution of the eigenvalue problem for the Hamiltonian ()3.27|) can be obtained within the 
method of the Baxter Q— operators 0. The standard method of solving spin chain models, the 
Algebraic Bethe Ansatz (ABA) method p3 EH] , is not applicable in the case under consideration 
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since the Hilbert space of the model does not have a normalizable lowest weight vector. 

We remind that by definition the Baxter Q— operator is the operator which acts on the Hilbert 
space of the model, depends on the spectral parameter and satisfies the following requirements 

• Commutativity 

[Q(u),Q(V)] = [Q(u),t N (v)] = [Q(u),H N ] = 0, (4.1) 

• Baxter equation 6 

t N (u)Q(u) = (u + is q ) N Q(u + i) + (u - is q ) N Q(u - i) . (4.2) 

In the case under consideration, the spectral parameter, u, as natural to expect, has the same 
nature as the spectral parameter of the TZ matrix, u = (u,e). We also used the notation i = 
(i, 1/2) and s g = 1/2 + ip g in Eq. (JOJ). 

As usual we look for the operator Q in the form of an integral operator. Then Eq. ()4.2jl results 
in a particular differential equation on the kernel of the operator Q. One can find the general 
solution of this equation using approach developed in [7] . This is based on the invariance of 
the auxiliary transfer matrix ()3.31|) with respect to local rotations of the Lax operators 

L k {u)^L k {u) = M^ 1 L k (u)M k+1 , (4.3) 

where M k are arbitrary 2x2 nondegenerate matrices. Then one can show that the function 

N 

Y u (x, y) = JJ D au (x k - y k+l ) Df3 u (x k - y k ) , (4.4) 
k=i 

where 

« u = (s g - iu, e + e q ) , (3 U = (s q + iu, e) , a u + (3 U = (2s q , e q ) , (4.5) 

satisfies the following equation 

t N (u)Y u (x,y) = (u + ts) N Y u+i (x,y) + (u-ts) N Y u ^(x,y) (4.6) 

at arbitrary values of variables yi, ■ ■ ■ ,Vn- This means that the convolution of the function 
Y u (x, y) with an arbitrary function of Z(y) also satisfies Eq. ()4.2|) . The details of the derivation 
of Eq. ()4.6|) can be found in Ref. ^l] where the similar case of the SX(2,C) spin chain was 
considered. Let us only note here that one can consider the function Y u (x, y) as an array of two 
functions which depends on the complex parameter u only, Y u = (Y( u ,o), ^(u,i/2))- Then Eq. f!4.6|) 
takes form of the matrix finite-difference equation. 

Thus the kernel of the Baxter operator can be written in the form 

f N 

^ fc=i 

6 Similarly to the Yang-Baxter equation 13.1J1 this equation contains two finite-difference equations which en- 
tangle the Baxter operators Q(u,0) and Q(u, 1/2). 



13 



Vi V2 y N Vi 2/2 y N 




Xi x 2 x 3 x N X\ Xi y x 2 y x 3 x N y X\ 



Figure 8: The kernel of the Q-operator in two equivalent representations. The indices are ex- 
plained in the text; x' = 1 — x. 

where the function Z(y, x') does not depend on the spectral parameter u. The restrictions on this 
function can be deduced from the requirement of the commutativity of the Baxter Q— operator 
and the auxiliary transfer matrix. To this end let us note that the function Y u (x, y) possesses 
the following property 

Y u (x,y) = (-l) 2m «VY u ,(y,x), (4.8) 

where u' = (—u, e + e q ) and V is the operator of the cyclic permutation ()3.35jl . Then taking into 
account Eq. ()4.6|) and the invariance of the transfer matrix under cyclic permutations, [P, tif(u)) = 
0, one derives 

t N (u,S(p g ,x))Y u (x,y) = (-l) N t N (-u, S(p q ,y))Y u (x, y) = t N (u, -S(p q ,y))Y u (x,y) . (4.9) 

Here we indicated explicitly that the transfer matrix tjv(w) is expressed in terms of the differential 
operators ()2.6|) acting on x or y-coordinates. Then, noticing that the integration by parts results 
in the change p q -> -p q in the generators / dy [S(p q , y)${y)] #{y) = J dy$(y) [S(-p q ,y)^(y)] 
and taking advantage of Eq. ([4.9)1 . one derives that the requirements [t N (u), Q(u)] = gives the 
following equation on the function Z(x, y) 

t N (u,S(-p q ,y)) Z(y,x") = Z(y, x') t N (u, S(p q , x 1 )) . (4.10) 

The simplest solution to this equation corresponds to the operator Z intertwining the generators 
Sk{—p q ) and Sk(pq)- Thus one can choose the operator Z to be proportional to the product 
M.% ® M.2 ® ■ ■ ■ ® M-n ', where the intertwining operators M.^ are defined in ()2.3jl . For later 
convenience we put 

N 

Z(x,y) = U^yM-Xfc). (4.11) 

fc=i 

with 7 = (2 — 2s q , e q ). 

As usual, we visualize the kernel of the Q— operator as a Feynman diagram, which is shown in 
the lhs of Fig.|H| In the rhs of Fig.|HJ we give the equivalent representation of the kernel ()4.7|) which 
can be obtained from the diagram on the lhs with the help of the star-triangle relation ([3. 10)1 . The 
graphical representation of the kernel is very convenient for the analysis of the properties of the 
Baxter Q— operator. As an example we prove the first relation in (j4.1|) . namely [Q(w), Q(v)] = 0. 

The diagrammatic representation of the integral kernel of the operator Q(u) Q(v) is shown 
in Fig. El There we represent the kernel of the operator Q(v) by the diagram shown in the 
lhs of Fig. |S1 while for the operator Q(tt) we choose the alternative representation given by the 
diagram in the rhs of Fig.|8j We also inserted two propagators with opposite indices in the central 
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Vi V2 Vn 




X\ X 2 £3 ^4 X N X\ x 2 

Figure 9: The diagrammatic representation of the product Q(w)Q(t> ). 




Figure 10: The diagrammatic representation of the exchange relation. The indices take the special 
values ol x = (s+ix, e x + e) and /3 X = (s — ix, e), where s = 1/2 + ip, (a u + (3 u = et v + f3 v = (2s, e)). 

rhombus in Fig. El Since Do-(x) D^^x) = 1, such an insertion does not alter the value of the 
diagram. 

The proof is based on the use of an additional diagrammatic identity- exchange relation- 
which is shown in Fig. El This identity can be easily derived with the help of the star-triangle 
relation ()3.10|) . Further, to prove the commutativity of the Baxter operators we choose the index 
of the inserted propagator to be, cr = (3 U — f3 v . Then one can use the exchange relation to move 
this inserted propagator, D^, to the right. The displacing of the propagator from left to right 
of the cross subdiagram alters the indices of the propagators and produces some scalar factor in 
the way shown in Fig. After repeating this operation N times the propagator resumes its 
initial position and annihilates the propagator D_ a . Therefore the resulting diagram will have 
the same form as that in Fig. El Examing the indices of the propagators and the scalar factor 
one finds that they correspond to the diagram for the operator Q(v ) Q(u), if the representation 
in the lhs of Fig. |H] is used for Q(n) and for operator Q(v) those in the rhs of the same figure. 
Thus one concludes that 

QWQW = Q(v)Q(u). (4.12) 

The Eqs. f!4.4|) . ()4.11|) together with f)4.7|) define the Baxter Q— operator which satisfies Eqs. (j4.1j) . 
()4.2|) . (Strictly speaking we did not show yet that the Baxter operator commutes with the 
Hamiltonian but will do it later on.) 
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4.1. Properties of Q— operator 

Let us examine the Baxter Q— operator for the special values of the spectral parameter such that 
indices a or f3 become equal zero. Using the diagrammatic representation of the latter given 
in Fig. |H1 one finds that for u = (is q , 0) (/3 = 0) the lines with the index /3 disappear and the 
integral over centers of the star diagrams can be easily calculated with help of Eq. ()3.9|) . A short 
calculation gives 

(4.13) 
(4.14) 



Q(is„0) = c(p q )"I, 
where I is the unit operator and the normalization constant is given by 



c(p) = n\A(l + 2ip,e) 



7T 



P 



tanh 46 - 1 ^) . 



Repeating this calculation for u = (—is q , e q ) (a. = 0), one finds 



c( Pq ) N V. 



. ■ (4.15) 

where V is the operator of the cyclic permutations ()3.35|) . Therefore, the operator V can be 
expressed as the ratio of the Baxter operator at two special points 



V 



"iS q , € q 



is q ,0) . 



(4.16) 



Next, using the diagrammatic representation for the Baxter Q— operator (Fig. |HJ) and for the 
transfer matrix (Fig. one finds after some algebra 



X(u, p , p g ) [Q{u* - is , e + eo)] 1 " Q{u + is , e) 

= x{u, Po> P g ) Q(u + i{l - s ), e + e ) [Q(u* - i(l - s ), e)]t 
The normalization factors are given by the following expressions 



X(u, p , p 



D iir(eo+e q )N 



-2iu 



71 



N 



x [A(l + ip q + ipo -iu,e + e q ) A(l - ip q - ip -iu,e + e )] 
(-l) 2 ^ )N c(p q )- N r PoPa (uf 



-N 



(4.17) 
(4.18) 

(4.19) 
(4.20) 



X( U , P0' Pq) — K^^J ~ ^KHqJ ' PoP q 

where r PoPq {u) is defined in Eq. f!3.16|) . Thus an arbitrary transfer matrix ()3.25|) can be expressed 
as the product of two Baxter Q— operators at special values of the spectral parameter. 

Using the diagrammatic representation, it is also straightforward to show that the operators 
Q and commute at arbitrary values of the spectral parameters. Moreover, one can deduce the 
following identity 

[(-1)*- A(a v ) A((3 V )] N Q(«) [QK)f = [(-l) 2eu A(a u )A(f3 u )] N Q(v) [Q(u*)]t , (4.21) 

where oc u , f3 u are defined in ()4.4|) . 

Having put p = p q in Eq. ()4.17|) and taking into account ()3.27|) . one obtains the following 
expression for the Hamiltonian in terms of the Baxter Q— operator 



H N = i^-\nQ(is q + u,0) 
au 



u=0 



+ i ^j-\nQ(-is q + u } e q ) ] 



+ £ 



N 



u=0 



.d 

% — In w 

du 



[is q + it, e q ) 



u=0 



+ i — In u 

du 



N a 



-is* + u, 0) 1 



u=0 



£ 



N 



(4.22) 
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Figure 11: Distribution of the poles of the Baxter Q-operator in the complex w-plane for 
e = e g = 0. 

where the additive constants £jy and £n are 

£ N = N [4j{2e q + 2ip q ) + ^{2e q - 2ip q )] , £ N = 2Ni/;(l). (4.23) 

Thus the Baxter Q— operator commutes with the Hamiltonian TIn (j3.27j) . Evidently, the rela- 
tions (|4.22jl . ()4.16|) . (jUTTj) hold for the eigenvalues of the corresponding operators as well. Thus, 
the knowledge of the eigenvalue of the Baxter Q— operator allows one to restore the eigenvalues 
of all other operators in question. The remarkable property of the Baxter Q— operator is that 
its eigenvalues can be determined without solving the eigenvalue problem. Indeed, since the 
Baxter Q— operator satisfies Eq. (|4.2)l . the same equation holds for its eigenvalues. Therefore 
solving this equation in the appropriate class of functions one can determine all eigenvalues of 
the Q— operator. In order to be an eigenvalue of the Baxter Q— operator the solution of Eq. ()4.2|) 
should satisfy the additional conditions. To find them we shall study in the next section the 
analytic properties of the eigenvalues of the Q— operator as the function of u. 

4.2. Analytic structure and asymptotic 

We shall denote the eigenvalue of the Baxter Q— operator corresponding to the eigenfunction 
^ q2 ,...,q N (%) as Qq(u) = Q q (w, e) and study its properties as a function of the complex variable u. 
It follows directly from the definition (j4.4j) and (j4.7J) that the Baxter Q— operator is well defined 
as the integral operator in the strip —1/2 < iu < 1/2. To continue it to the whole complex 
plane, let us apply the Baxter operator, Q(u), to a test function Then it can be seen 

that the resulting function = [Q(w)^](x) as a function of u admits continuation to the 

whole complex plane except for the points where the indices a and (3 take positive integer values, 
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Lower half-plane 



namely u~ = is q — in, (/3 = n > 0) and u+ = —is q + in (a = n > 0). The structure of the 
singularities can be easily established using the representation of the Baxter Q— operator given 
in the rhs of Fig. |H| Since the indices of the propagators a' and f3' are not positive integers 
while a or (3 are, the corresponding integrals define a regular function of u at these points and 
all singularities are contained in the prefactor A(ct, f3, *y) N ■ The straightforward analysis shows 
that the Baxter Q— operator has poles of the order N in the upper and lower half-planes in the 
following points 

• Upper half-plane 

e = eg, ul = p q + 2m + I, n>0 

(4.24) 

e ^ E q , ul = p q + 2m, n > 1 

e = 0, u~ — —pg — 2in — |, n > 

(4.25) 

e = 1/2, u n =-p q -2in, n>± 

Thus the eigenvalue of the Baxter Q-operator is a meromorphic function of the variable u with 
poles of the order iV at the points «^(e), Eqs. (|4.24|) and ()4.25|) . The positions of the poles 
in the u plane are schematically shown in Fig. ^2 Note also that the equation in the second 
line of ()4.22j) gives the energy as the ratio of the residues at the leading (of order N) and the 
subleading (of order N — 1) poles at = ±(p q + i/2). 

Since the solutions of the Baxter equation (|4.2|) are determined up to multiplication by the pe- 
riodic function f(u+i) = f(u) one needs additional information to fix this arbitrariness. This can 
be obtained by the study of the asymptotic behaviour of the eigenvalues at large u. First of all, let 
us notice that the transformation properties of the eigenfunction under SL(2, M) rotation is deter- 
mined by the eigenvalue of the charge q2 or the total Casimir operator, 5* 2 = ^Si + . . . + S 

(see Eq. ()3.33J) ). Namely, S 2 = — 1/4 — p 2 if the eigenfunction transforms according the irre- 
ducible representation of the principal series, T^ p,e \ and S 2 = h(h — 1) if it transforms according 
to the discrete series representation V^. Thus depending on the value of the Casimir operator, 
the eigenfunction can be represented as follows 



2 

N 



%( Xl ,...,x N ) = / dx U q P± qN (x 1 ,...,x N \x )f ip ' e \xo), (4.26) 

J — oo 

^( Xl ,...,x N ) = Jv ± wn^g N (x 1 ,...,x N \w)<t>±(w), (4.27) 

where the integration measure in (jPZD is defined in Eq. (TP1) . The projectors D>> e ), 
intertwine the representations (<&T Pq ) N and T^ p,t \ T>^, respectively. They have the following 
transformation properties 

k=l 
k=l 
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where x' k = (ax k + b)/ (cxk + d) and w' = (aw + b) / (cw + d)\ a, b, c, d are real and ad — be = 1. 
In particular, one finds that the projectors are invariant under the simultaneous shift of all 
arguments by a real number and are transformed under scale transformations as 

n^ )to (Ax 1 ,...,Ax JV |Axo) = \- 1+s - Ns m^l qN (x 1 ,...,x N \x ), (4.30) 
Xx N \\w) = A * V -'I ij;; , /x (,r ; v N \w) . (4.31) 

Applying the Q— operator in the form (|4.7jl to the eigenfunction ()4.26|) . one finds that the leading 
contributions at u — > oo comes from two integration regions over y 

(I): \y k \=0(u), (II): y fc -2/ fc+ i = 0(l/«). (4.32) 

Next, let us notice that any function f^ p,t \x) can be represented as the transformation of the 
function fo(x) = 1, p p ' e '(x) = J Dg(p(g)T^ p ' e \g)f (x), where the integral is taken over the group. 
Then taking into account the invariance of the Baxter Q— operator with respect to the SL(2, M) 
transformations and the properties of the projector H^fl, qtf (|4.28j) . one can derive (see Ref. [TT] 
for details) 

Q q (u,e) u =°° (A lU s - Ns « + An u 1 ~ s ~ Nsq ^ [1 + 0(l/u)) . (4.33) 
The constants Aj and An depend on the integrals of motion q^, and we remind that s = 1/2+ip, 

S q = 1/2 + p q . 

In the case of the eigenf unctions of the discrete spectrum (|4.27J) . a careful analysis shows that 
the contributions coming from the regions (/) and (II) (|4.32|) are of the same order 

Q q (u,e) u =°° Cu 1 - 11 -^ [1 + 0(l/u)] . (4.34) 

Let us also note that since the kernels of the operators Q(u, 0) and Q(u, 1/2) coincide in the 
regions (J) and (II), the difference Qq(u, 0) — Qq(u, 1/2) vanishes faster than any degree of 1/u, 
i.e. 

Q q (u, 0) - Q q (u, 1/2) ^=°° 0(e-^) , (4.35) 

where k is some constant. 

Thus the solution of the Baxter equation ()4.2j) corresponds to some eigenvalue of the Baxter 
Q-operator only if it has the proper pole structure ()4.24|) . ()4.25|) and proper asymptotic (|4.33)L 
(PDMl) at u — ► ±oo. 



5. Separation of Variables 

In this section we explicitly construct the representation of the Separated Variables for the model 
in question. Namely, we shall obtain the following integral representation for the eigenfunctions 
of the model 

*q(^) = 0~ dp JJ^fcAtJv(^) < ^q(P ) ^l,---,^7V-l)f / P, 1 xi,...,'Uiv_ 1 (^), (5.1) 
71 J k=l 

where 

S Vu *=w?jZ dUk 5, ' (5 - 2) 



£fc=0,l/2 
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/xjv(w) is the Sklyanin measure (see Eq. (j5.12jl ). $ q (p, «i, • • • , Mjv-i) is the eigenfunction in 
the SoV representation (j5.16|) and U PtUlt _,^ UN _ 1 (xi, . . . ,xn) is the transition kernel to the SoV 
representation (J5.5J) . 

In Sklyanin's approach jl] the functions Up^^^^xi, . . . ,xjy) are identified as the eigen- 
functions of the operator B^{u) — an off diagonal matrix element of the monodromy matrix. 
Provided that the spectrum of the operator Bn{u) is non-degenerate, one can derive that the 
eigenfunction in the SoV representation ($ q (p, tti, . • . , Uzv-i)) satisfies the Baxter equation in 
each variable. In the case under consideration this condition is not fulfilled. Moreover, the 
hermitian operator B^iu) admits different non-equivalent self-adjoint extensions and one has to 
choose the correct one. Although these problems can be overcome, we shall construct the tran- 
sition kernel U PiUlt ,,, tUN _ 1 (xi, . . . ,xn) following another approach. It was conjectured in Ref. j22] 
that the transition kernel to the SoV representation can be related to the kernel corresponding 
to the product of the Baxter Q— operators. The representation of such type is known now for 
a number of spin chain models fTJ ■ We suggest the following ansatz for the transition 

kernel 

N 

U P:Uu ...^ N _ 1 {x) = cjv(p) Yldy k e ipm [Q(ui) . . . Q(wjv-i)] fa, ■ ■ ■ , x N \yi, . . . , y N ) (5.3) 

fc=i 

and show that so constructed kernel possesses all necessary properties. For convenience we put 

c N {p) = \p\~ N+l/2 c{p q )- N ^l\ 

where c(p q ) is defined in (I4.14j) . 

First of all let us notice that, as follows from the properties of the Baxter Q-operator, the 
kernel ()5.3|) is a symmetric function of the variables Ui, . . . , u N _i. Second, it has a definite mo- 
mentum i(S_ + • ■ • + S^)U Pi a(x) = pU P} a(x) and, at last, it satisfies the Baxter equation ()4.2|) 
in each variable u^. 

Further, we shall show that for real k = 1, . . . , N — 1, the functions U p ^(x) are mutually 
orthogonal 

J l[ dx k (U g ,t(x))* U p ,4x) = 2n5(p - q) 5{u - v) , (5.4) 

where 

5(U - v) = ^ 5(ui - v ix ) . . . 5(wiv_i - v iN _ x ) , 
s 

the sum goes over all permutations and 5{u — v) = (27i) 2 5 tuev 5(u — v). To prove ()5.4|) it is 
convenient to represent the kernel (J5.3)) in another form. Namely, let us use the diagrammatic 
representation for the Baxter Q-operator shown in the rhs of Fig.lHl The diagram for the product 
of N — 1 Q-operators consists of N — 1 such rows. We shall assume that each triangle in Fig. |8] 
accompanied by the corresponding scalar factor ^A(cx, f3, 7). According to Eq. ()5.3|) one has 
to perform the integration over variables y2, ■ ■ ■ ,Vn which correspond to the upper vertices of 
the N — 1 triangles in the upper row. Since only two propagators are attached to each vertex 
in question, one can use the chain relation, Eq. ()3.9|) . The resulting propagators have the index 
7 = (1 — 2s q , e q ) and cancel the N — 1 horizontal lines in the upper row of the diagram. Thus 
after the integration over its upper vertex, the triangle disappears producing the scalar factor 
c(p q ) = ir\A((l + 2ip q ,e q ))\ 2 . As a consequence one obtains a new diagram in which N — 1 
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Figure 12: The graphical representation of the transition kernel, Eq. ([5.5)1 . 



triangles in the upper row are eliminated and the additional scalar factor c(p (? ) iV_1 is acquired. 
Repeating this operation N — 1 times one derives that the kernel ()5.3|) is represented by the 
diagram shown in Fig. ^] with the accompanying additional factor c(p 9 ) 7V( ^ iV_1 ^ 2 . Therefore the 
function ()5.3|) can be represented in the form 

= \p\- N+1/2 J dye*™ [A N (u 1 )...A 2 (u N _ 1 )](x lj ... J x N \y), (5.5) 

where the function A k (xi, . . . , Xk\yi, ■ ■ ■ , yk-i){u^ + i^ k ) corresponds to the k — 1— th row (from 
the top) of triangles in the diagram shown in Fig. We remind here that each triangle is 
accompanied by the corresponding scalar factor. 

The representation ()5.5|) is convenient for the calculation of the scalar product ()5.4|) . Using 
the representation ()5.5|) for both U— functions, one finds that the scalar product ()5.4|) takes the 
form (up to some prefactor) 

dy' e**t j dy e™ [A 2 { Vl )1 . . . A N (v N ^y A N { Ul ) . . . A 2 (u N ^)] (y'\ y) . (5.6) 

It can be easily checked that for real u, v, (u ^ v), the A— operators satisfy the following 
(exchange) relation 

A k (vyA k (u) = c(p g )<p(u,v)A k _ 1 (u)A k _ 1 (v)i (5.7) 

with 

n t_\\2e u +2e v 

^ u ' v) = A(0,-0,)A( a .- a .y (5 - 8) 
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Using this relation one obtains for ()5.6|) 

N-2N-l-j 

c(p q )<»-w-w n n vfavH-j) (5.9) 

3=1 k=l 

x J dy'e-^' J dye*™ [A 2 (v x y A 2 (u 1 ) . . . A 2 (v N ^A 2 (u N ^)] (y'\y) , 

where we suppose that u k 7^ Vj if k ^ j. The calculation of the convolution of the A functions 
in the second line of ()5.9|) (which can be represented as a chain of the box diagrams) is relied on 
the identity 

ipx 

dy^[A 2 {vfA 2 {u)]{x\y) = c (p q )— (2tt) 2 6 w J(u - v) . (5.10) 
Collecting all factors together, one derives for (15. 6|) 

N-2 N-l 



c (Pc 



vJV(AT-l)/2 



11 II V{UrmUj) 



m=l j=m+l 



N-l 



(2tt)5(p - q) Y[ S(u k - v k ) 



(5.11) 



k=l 



where 8(u — v) = (2ii) 2 5 euev 5(u — v). Restoring the symmetry with respect to the permutations 
of variables u±, . . . , Un-i (vi, . . . , Ujv_i), which is broken due to the imposed condition, u k 7^ Vj 
for k 7^ j, one arrives to the expression in the rhs of Eq. ()5.4|) with 



(N-l)\ 



-N(N- 



N-2 N-l 

1,2 n n 

m=l j=m+l 



UJ{U m ,U jy 



where 



-1 1 l^j 



tanh7T 1 



l-4(e m +ej) 



(5.12) 



(5.13) 
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We remind here that (e m + e,) is equal to zero if e m = €j and to 1/2 otherwise. 

The measure (1n(u) is a regular nonnegative function for real u k . When one of the separated 
variables u k goes to infinity, it grows as 



/ -,\ Mfc^dzOO I iN-2 



i + 0(i/M)) 



(5.14) 



Under the analytic continuation to the complex plane, the measure ()5.12j) becomes a meromorphic 
function of the variables u k . One can also verify that it satisfies the functional relation 



flN{Ul, ...,U k + l,... ,UN-l) 
HN(u h . . ., life, . . . ,ltjv-l) 



n 



u k —Uj + i 

U k - Uj 



(5.15) 



To determine the wave function in the SoV representation, ^ q2 ,...,q N (p, u%, . . . , it/v-i), one cal- 
culates the scalar product of the eigenfunction in the coordinate representation, tyq(x), (Eq. 1 )5. 1^ 
with U„ 



'p ) t*i,...,«jv-i(^)> ( E 1- (ESI)- Using Eq. (}5~4")) one finds 

$ q (p, Ul, . . . ,U N -l) = Cy{p) Qq(Wl)* . . . Qq^jV-l)* 



where 



c*(p) = catQo) / dx 1 ...dx N e ipxi ^q(xi, . . . , xtv) . 



(5.16) 
(5.17) 
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Assuming that the basis formed by the functions t/ p ,u 1 ,..., Ujv _ 1 (x) is complete, one derives the 
orthogonality condition for the wave functions in the SoV representations 

/N-l N-l 
Y[ T>U k fi N (u) JJ Qq'(Wfc)Qq(Wfc)* ~ 5qq' , (5-18) 

k=l k=l 

where we took into account that the integral over p is factorized. 



6. Special case: N = 2 



In the case of the two-sites spin chain, the eigenvalues of the Baxter Q— operator can be obtained 
in a closed form. Indeed, for the N = 2 spin chain, the eigenf unctions are fixed by the group 
properties, Eqs. ([A.4J1 . (|A.9[) . Henceforth, for brevity, we shall restrict ourselves to the case of 
the representation of the positive parity, i.e. we shall imply that e q = 0. Applying the Baxter 
Q— operator to the eigenfunctions and going through the calculations, one finds the represen- 
tation of the Baxter function in the form of a one-dimensional integral. We shall denote the 
eigenvalue of the Baxter operator corresponding to the eigenfunction of the continuous series 
(n^' e p (x, yi, 2/2))*) (see Eq. (IA.4jl ) by Q e p (u, e) 7 . The eigenvalues of the Baxter operator corre- 
sponding to the eigenfunctions of the discrete series are the same, Q(u, e)(U p ^ p ^(w, xi, x 2 ))* = 
Qh(u, e)(n.p'p 2 (w, x\, X2))* . Up to unessential here u— independent factors the eigenvalues of the 



Baxter operator are given by the following expressions 

Qh(u, e) = Q h (u,e) + (-l) h Q h (-u,e) 



where 



Q%u,e) = Qt p ( u ,e) + (-iyW p (-u,e), 

Qh(u, e) = / drr lu - 1 [q h (r) + (-l) 2e q h (-r 
Jo 

Q%u,e) = f dTT™- l [q%T) + {-iy%{-T)] . 
Jo 



(6.1) 
(6.2) 

(6.3) 
(6.4) 



The functions qh and q £ p can be expressed in terms of the Legendre functions of the second kind 
as follows 



<f {r) 



'1-T 



> 1 Hh-l 



1 - T 



1- 



\{p,e) Q_i 



where 



X(p,e) 



/2+ip 



l + (-l 



1+r 



,2e 



A(-p,e)Q. 



1/2-ip 



T 



T 



sinh up 



(6.5) 



, (6.6) 



(6.7) 



7 Since we assume that p q = (p q ,0), the irreducible representation of the positive parity (e = 0) appears in 
the tensor product decomposition, therefore p = (p, 0). We remind also that the parameter e = 0, 1/2 marks two 
different eigenfunctions having the same value of the two-particle Casimir operator. 
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The argument of the Legendre function entering the integrals 1)6.3)1 . (J6.4j) varies in the range 
[— 1, oo]. We remind here that for real x, — 1 < x < 1, the Legendre function is defined by the 
relation Q„(x) = |(Q I/ (x + it) + (Q v (x — ie)). 

The factors (— l) h , (— l) 2e entering Eqs. ()6.1)1 . (J6.2)) determine the parity of the eigenfunctions 
with respect to the permutation of the arguments. Therefore, one concludes that the eigenvalue 
of the Baxter operator is even (odd) function of u for the eigenfunction of positive (negative) 
parity. 

Further, it is evident that the function Q{u) (Q(—u)) has poles in the upper (lower) half-plane. 
The poles of Q{u) occur at the points where the integrals f|6.3|) start to diverge at r — > 0. Noticing 
that the integrand in (J6.3|) . (J6.4|) is an even (e = 0) or odd (e = 1/2) function of r and taking 
into account the properties of Legendre functions, one concludes that the function Q(u,e) has 
poles of the second order located at the points ()4.24jl . 1)4.25)1 . It is also straightforward to see that 
the asymptotic of the Baxter functions 1)6.1)1 . ()6.2)1 at u — > oo is in agreement with Eqs. ()4.33j) . 
f!4.34j) . Calculating the difference Q(u, 0) — Q(u, 1/2) one finds that it can be represented as 

POO 

~ / drr^qi-r) , (6.8) 
Jo 

where q{r) is given by (J6.5j) or ()6.6)) . Since the function q(r) is analytic near r = — 1 we reproduce 

At last, to find the energies of the corresponding eigenstate it is necessary to calculate the 
expansions of the Baxter functions near the first poles, u = ±i(l — s q ) (see. Eq. (J4.22)) ). Taking 
into account that for r — > 

(1^7) = ~\ l0g |r| + - ^(1 + ^) + 0(t) (6.9) 

one finds, e.g. for Q £ (—i{l — s q ) + u, 0), 

- f 1 dr 

Q e p (-i(l-s q )+u,0) = 2 - 1 _[-logr+ (6.10) 

^0 T 

A(p, e) [^(1) - ^(1/2 + z»] + A(-p, e) [^(1) - ^(1/2 - ip)]] + . 

Then using Eq. ()4.22)) and taking into account that 7i 2 = 2if 12 , one reproduces the expres- 
sions icragi) . 

7. Summary 

In this paper we studied the spin chain model with the SX(2,R) symmetry group. The Hilbert 
space attached to each site is L 2 (M) and the symmetry transformations are realized by the 
operators of the unitary principal series continuous representation of the SL(2, R) group. To 
define the model we constructed the SL(2, R) invariant solution of the Yang-Baxter equation 
- the TZ— operator which acts on the tensor product of two L 2 (R) spaces. The Hamiltonian of 
the model is defined as the derivative of the fundamental transfer matrix and is given by the 
sum of the two-particle Hamiltonians. The pair-wise Hamiltonians have both the discrete and 
continuous spectrum, Eq. (|3.29jh which reflects the pattern of the decomposition of the tensor 
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product of two representations of the principal continuous series into irreducible components. 
The eigenstates belonging to the continuous spectrum are specified by the value of the si (2) spin, 
s = 1/2 + 2/?, and parity with respect to the permutation of the coordinates x\ and x-i. The energy 
gap between the eigenstates with the same value of the spin and different parity is maximal for 
p = and exponentially decreases with p. 

To solve the model we applied the method of the Baxter Q— operator p] and the Separation 
of Variables |IJ . The standard ABA approach [2] is not applicable for the model in question due 
to the absence of the lowest weight vector in the Hilbert space of the model. Having realized 
the Baxter operator as an integral operator we resolved the defining equations and obtained the 
kernel in the explicit form. It allowed to determine the properties of the eigenvalues of the Baxter 
operator as functions of the spectral parameter. Then the eigenvalues of the Baxter operator 
can be obtained as the solutions of the Baxter equation in the certain class of functions. The 
eigenvalue of the Baxter operator encodes all information about the corresponding eigenstate. 
We have shown that the Hamiltonian of the model can be obtained as a derivative of the Baxter 
operator at special points. Moreover, the arbitrary transfer matrix factorizes into the product 
of two Baxter Q operators at certain values of the spectral parameters. Analoguous results have 
been obtained for the noncompact SX(2,R) (discrete series) and SL(2,C) (continuous series) 
spin chains, see Refs. [TT] . 

We have constructed the representation of the separated variables for the model in question. 
The kernel of the unitary operator, which maps the eigenfunction to the SoV representation, has 
been obtained in an explicit form. It factorizes into the product of iV — 1 (N is the number 
of sites) operators each depending on one separated variable only. The kernel of the transition 
operator can be visualized as a Feynman diagram with a specific pyramidal form. This form of 
the kernel, first obtained for the SL(2, C) spin chain ill , is a general feature of all noncompact 
SL{2) spin magnets [121 US]. Using the diagram technique we calculated the scalar product of 
the transition kernels and determined the Sklyanin's integration measure. We have shown that 
the wavefunction in the separated variables is given by the product of the eigenvalues of the 
(conjugated) Baxter operator. Therefore the knowledge of the eigenvalue of the Baxter operator 
allows to restore the eigenfunction. 
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A Appendix: Tensor product decomposition 

In this Appendix we collect the necessary formulae concerning the decomposition of the tensor 
product of the representations of the principal continuous series. The tensor product of two 
representations of the principal continuous series is decomposed into the direct integral over the 
representations of the same type and into the direct sum over the unitary representations of the 
discrete series (2^]- We remind that the representations of the discrete series T>^ are labelled by 
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the (half)integer number, h, and can be realized by the unitary operators D£(g) (D h (g)) 

[Dt(9)*±] H = {cw + d)-»9 ± (^) , (A.l) 

acting on the Hilbert space H^, respectively. The latter are defined as the space of the functions 
analytic in upper (TC^), lower (TC^) half-plane with the scalar product defined by the following 
expression j25] 



(*±|$±) = J V±wV±(w)®±(w), (A.2) 

where 

oh _ i 

V±w = e(±lm(w))(±lm(w)) 2h - 2 dxdy (A.3) 

7T 

and w = x + iy. 

The operators separating the irreducible components in the tensor product T Pl ® T P2 are 
the following: the projectors to the discrete series, n^^, which map TPi TP* -> D±, and the 
projectors to the principal continuous series. The projecting operators can be realized as 

integral operators. The irreducible representation T p , (p = (p, e), p > 0, e = e\ + e 2 ) enters 
into the tensor product T Pl (g) T P2 with double multiplicity. Therefore one can construct two 
projectors, Rpf P2 , where the parameter e takes two values and 1/2 and marks these equivalent 
representations. The integral kernel of the projector is given by the product of three propagators 

n pf P2 (3,2/1,2/2) = D a3 (y 2 -x)D a2 (x -y 1 )D ai (y 1 -y 2 ) , (A.4) 
where the indices are defined as follows 

<*i = (g ~i(P + Pi + &),£ + d), 

«2 = d +i(p + P2 -pi),e), (A.5) 
c*s = (~ + «(p + Pi -p 2 ),£ + ei + £2) • 
Then the function $ p,e defined as 

$ p ' £ (x) = J dy 1 dy 2 W 3 pi £ p2 (x,y 1 ,y 2 )y(y 1 ,y 2 ), (A.6) 

transforms according to the representation T p . Next, using the techniques for the calculation of 
the SL(2, R)-integrals introduced in Sec.EU it is straightforward to check the following orthogo- 
nality condition 

J dy 1 dy 2 n% p2 (x,y l ,y 2 ) (u^{x' ,y u y 2 ))* = N{p)5 ee , 8(x - x')8(p- p'), (A.7) 

where 

M{p) = (27r)V 1 coth 1 - 4e ( 7 rp). (A.8) 
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The conformal spin h of the representations of the discrete series entering the decompo- 
sition of the tensor product T Pl ® T P2 is integer for e\ = e 2 , and half-integer otherwise. The 
projecting operators onto irreducible components are 



Ilh pfp 2 {w,x 1 ,x 2 ) 



(x 2 - xi) 



h+ei+E2 



X 



X\ — W / 

X\ — w 
X 2 — w 



h+e 2 -eif X2 _ w 



(Xi - X 2 ) 



h+ei-t-2 
[i(pl-p2)~ei+e2]/2 



Xx - X 2 



(A.9) 



x 2 — w 
X\ — w 



(x 2 - Xi) 



[i(/?2— pi)— 62+ei]/2 



where the w lies in the upper half-plane for the "+" projector, and in the lower for "— " projectors. 
Note that the expressions in the square brackets in (jA.9|) are the single-valued functions of w in 
the upper or lower half-plane. Thus the functions $ h,± 



^(x) 



d yi dy 2 nj± 2 (w, y u y 2 ) y 2 ) (A.10) 
transform according to the representations T>^. The normalization condition reads 

J dx 1 dx 2 (ll^ P2 (w,y h y 2 )^ H^fay^yz) = c(h) 8 hN K(z, w), (A.ll) 
where the w) = e t7Th (z — w)~ 2h is the reproducing kernel (unit operator on Ti.^) and 



c(h) = (2tt 



i2 T(2h-l) 



T 2 (h) 



(A.12) 



Then any function ty(xi,x 2 ) G L 2 (M) (g) L 2 (IR) can be decomposed as 

00 „ 

*(ar 1 ,x 2 ) = J2 c-\h)jV ± w(n h pi ± p2 (w,y 1 ,y 2 )y^ ± (w) + 



h=l + (e 1 +e 2 )/2 
'■00 



dpM-\p) / dx (ll^ix^ux,)) $" £ (x), (A.13) 

where the functions <J) /l ' ± and <£> p,£ are defined in Eqs. (jA.10|) . ()A.6|) . respectively. 

Note, that if p l = p 2 , the projectors have definite parity with respect to the permutation of 
the coordinates X\ and x 2 , 



and 



. 1)/l+2ei n j± ( w ,x 2 ,x 1 ). 
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